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Two-dimensional AKLT model on a honeycomb lattice has been shown to be a universal resource for 
quantum computation. In this valence bond solid, however, the spin interactions involve higher powers 
of the Heisenberg coupling {Si ■ Sj)^, making these states seemingly unrealistic on bipartite lattices, 
where one expects a simple antiferromagnetic order. We show that those interactions can be generated 
by orbital physics in multiorbital Mott insulators. We focus on t 2 g electrons on the honeycomb lattice 
and propose a physical realization of the spin-3/2 AKLT state. We find a phase transition from the 
AKLT to the Neel state on increasing Hund’s rule coupling, which is confirmed by density matrix 
renormalization group (DMRG) simulations. An experimental signature of the AKLT state consists of 
protected, free spins-1/2 on lattice vacancies, which may be detected in the spin susceptibility. 

PACS numbers: 75.10.Nr, 75.10.Jm, 75.10.Hk 


Introduction - Spin liquids are ground states of mag¬ 
netic systems which do not exhibit magnetic symmetry 
breaking down to zero temperature [1]; in contrast to 
conventional ordered states like ferromagnets and antifer- 
romagnets (AF), they host fractional excitations [2]. A 
classical example is the spin-1 Heisenberg chain which, as 
predicted by Haldane [3, 4] and observed in experiment 
[5], does not order and exhibits fractionalized spin-1/2 
excitations at its boundaries. 

In 1987, Affleck, Kennedy, Lieb, and Tasaki (AKLT) 
formulated a family of exactly solvable valence-bond-solid 
(VBS) models which (i) account for the Haldane phase 
in one dimension [6], but also (ii) generalize to higher 
dimensions [7]. In these VBS states, the spin-S' at each 
site is partitioned into 2S spins-1/2, and the latter form a 
solid of singlet valence bonds, see Fig Id for the S = 3/2 
case. As a result the spin-singlet ground state exhibits 
spin-1/2 excitations at edges or boundaries, whenever 
valence bonds are broken. The model Hamiltonians which 
realize VBSs as exact ground states consist of a sum of 
operators projecting the total spin on neighbouring sites 
to the maximal possible value of 2S. These projectors 
may be expanded in powers of {Si ■ Sj); for the AKLT 
model of spins-3/2 on the honeycomb lattice which will be 
considered here it reads: 

Haklt = • Sgf + C3{S, ■ S,)% (1) 

(Lj) 

with C2 = HI and C3 = HI- The multi-quadratic in¬ 
teractions are typically not found in real materials e.g. 
insulators; usually the dominant spin-spin interaction is 
generated by superexchange, leading to AF Heisenberg 
interactions. Consequently, except for the ID case in which 



Figure 1. (Color online) Nearest-neighbout t 2 g hoppings, cor¬ 
responding to (a) xy, (b) yz, and (c): zx orbitals, and (d) the 
resulting directional-hopping model, (e) free spins (denoted by 
stars) around lattice vacancy in the AKLT phase. 

the absence of magnetic ordering at T=0 is guaranteed by 
strong quantum fluctuations, the ground states on bipartite 
lattices break the spin SU{2) symmetry and order antifer- 
romagnetically [7]. As a result, a physical realization of 
the 2D AKLT phase in bipartite systems remained elusive. 
The importance of just such states has been underlined 
recently by the discovery that they are a universal resource 
for measurement-based quantum computation [8]. 

In this paper we study theoretically multiorbital insula¬ 
tors described by Hubbard models with orbital-dependent 
directional nearest-neighbour hopping, and propose that 
they are a platform for realization of two dimensional VBS 
states. The choice of honeycomb geometry is motivated by 
recent activity on the layered hexagonal Iridates such as 
Na2lr03 and Li2lr03. In these spin-orbit entangled Mott 
insulators the t 2 g electrons lead to [9] anisotropic spin-spin 
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interactions of the Kitaev-Heisenberg model. We consider 
a similar Hubbard model, however in a d^ configuration, 
i.e. having three electrons per site (half filling); an example 
is Li 2 Mn 03 . While Li 2 Mn 03 shows magnetic order at low 
temperatures [ 10 ], variants of this material may be in the 
VBS phase. 

We show that the VBS state is realized in our model 
and is stable against interactions; it is adiabatically con¬ 
nected to the AKLT state. Using analytic arguments and 
DMRG simulations [11] we find that while strong Bund’s 
rule coupling leads to bulk Neel order, an extended region 
in the phase diagram at small Jh exists, with properties 
identical to the AKLT state (see Fig. 2): an energy gap, 
exponentially decaying correlations, a singlet groundstate, 
and, notably, symmetry protected free spin- 1/2 excita¬ 
tions at lattice vacancies, which provide an experimental 
signature. 

Model - We consider a tight-binding model of t 2 g elec¬ 
trons on a honeycomb lattice, arising for example in the 
basal plane of the Iridate compounds, as seen in Fig. 1: 

H = —t 

^ — T I'l' T 

+E E (2) 

i a' i 

where 7 = xy,yz,zx = 1,2,3 is the orbital index, 
^77 — ^ 2 , 7 ,(t^ 77 ,^ and Si — ^77 with — 

X^CTcr ''"17 (T^'^^G. 7 ,cr'- The notation implies that 

the hopping between nearest-neighbours i and j involves 
7 orbital only. 

The directional hopping originates from the anisotropy 
of the t 2 g orbitals, which are predominantly confined to the 
xy, yz, and zx planes. Consequently there is an appreciable 
spatial overlap of wavefunctions on neighbouring sites along 
preferred directions only, as shown in Figs, l(abc), and 
the hopping matrix elements in nonequivalent directions 
on the honeycomb lattice are dominated by distinct t 2 g 
orbitals, see Fig. 1(d). For simplicity we neglect hopping 
elements other than those in Eq. (2), which we don’t 
expect to change the picture. 

The interaction terms contain a Bund’s rule coupling 
favoring large spin formation, and a Hubbard interaction 
U.yY = U6.yY + U'{1 — S.yY). Herc U and U' are the 
interactions between electrons residing in the same or 
distinct orbitals, respectively; we assume U' < U. We 
focus on a configuration with 3 electrons per site, as in 
Li 2 Mn 03 [ 10 ]. 

Phase diagram - Consider first the case of vanishing 
Bund’s rule coupling Ju = 0. In the noninteracting limit 
U = U' = 0 the system trivially forms a dimer band- 
insulator (six electrons per unit cell); the ground state is 
a product of spin-singlets on links with an energy gap of 
2t. Allowing U ^ 0 while keeping U' = 0 does not alter 
the picture, since the Hamiltonian only couples pairs of 


Ujt 



Figure 2. Phase diagram of the system Eq.(2) for U' < U. 

sites connected by a single orbital 7 , i.e. it factorizes into 
a sum of decoupled dimer Hamiltonians One can 

show that the ground state remains a product of singlets 
for any value of U/t. 

In the large U limit three spin-1/2 states are occupied at 

each site, in distinct orbitals, since U' < U. The electron 

hopping generates an interaction between the spin degrees 

of freedom connected by the same orbital 7 , as shown in 

Fig. 1(d): Heff = JT,{i,j)^ Si .7 • S 77 with J = |^. Thus 

at Jh = 0 , for any {/, the ground state is a product of 

singlets with an energy gap Egap = 2\t\ ior t ^ U and 
,2 

Egap = jj = J for U t. While one cannot solve the 
model analytically for intermediate U, we expect a gapped 
phase all along the line Jr = 0, as illustrated in Fig. 2. 

Furthermore, at large Ujt and Jh U we can write 
down an effective Hamiltonian for three spins-1/2 per site, 
as a function of Jr'- 

Heff = J ^ ^ Si^ ■ Sjj — Jh 'y ^ ( ^ ^ 5 * 7 ) . (3) 

(Li>7 3 7 

In the limit of large Jn/t a, spin-3/2 degree of freedom 

is formed on each site i, denoted by Ti = Si,y They 

are coupled via a Heisenberg interaction H = JesTi ■ Tj 

with Jeff = f = Ifj) since in the spin-3/2 subspace of the 

Hilbert space of three spins-1/2 we have Sty = ^Ti. The 

spins-3/2 live on a bipartite lattice and will thus order 

antiferromagnetically in this limit and break the spin SU(2) 

symmetry. Since at small Jn/t the system is a gapped 

singlet, there is necessarily a phase transition as function 

■ 2 

of Jh- The transition should occur at Jr ^ J = jj- 
The spin-3/2 formation holds at large Jr for any U. In 
the limit oiU Jr their effective coupling is Jeff = 
(from 2 "^^ order perturbation theory). At 17 = 0 the full 
model ( 2 ) contains only two parameters {JH,t), hence we 
expect that the phase transition will occur at Jr ~ t. A 
schematic phase diagram is drawn in Fig. 2. We support 
it using DMRC simulations. 

Numerical simulations - To study the phase diagram 
we performed density matrix renormalization group (hnite 
DMRC) calculations using the ITensor (http://itensor.org) 
package. We simulated model Eq. (3) on a decorated 
honeycomb lattice using cylindrical boundary conditions 
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Figure 3. (Color online) (a) Energy gap of model Eq. (3) in 
units of J with and without an impurity. Inset: Finite size 
scaling of Egap for Jh/J ~ 10, (b) Inverse correlation length 
extracted from the exponential decay of (s^) (see Fig. 4a). 


Figure 4. (a) Expectation value on one of the sublat¬ 

tices, along the cylinder, 8 hexagons in length (b) Entanglement 
entropy Se of the ground state of model Eq. (3) normalized by 
the number of cut bonds. 


and system sizes of up to 216 spins-1/2, keeping a fixed 
aspect ratio of the cylinder. 

The calculated energy gap is shown in Fig. 3a: the 

system is gapped not only in the trivial case Jh = 0, where 
,2 

the gap is equal to ^ = J, but also over a finite parameter 
range of Jh < J- Upon increasing Jh the system reaches a 
transition at Jh J■ For finite size systems we simulated, 
the energy gap saturates at a finite value in the Neel phase, 
corresponding to the finite size spectrum of spin wave 
excitations. It vanishes, however, in the thermodynamic 
limit, consistent with the infinite size extrapolation shown 
in the inset. At small Jh the gap behaves as Egap = 
J — Jh + 0{Jjj). This linear dependence can be derived 
analytically; consider an elementary excitation of the VBS 
state, which is a triplet on some link with an energy cost 
Egap = J- In 1** order perturbation theory Jh induces 
hopping of such excitations on the Kagome lattice formed 
by the links, hence reducing the energy gap as shown. The 
phase transition may be thus thought of as condensation 
of such spin-triplet excitations. 

While in principle it is not obvious that there is a direct 
AKLT-Neel transition, the absence of an intermediate 
phase can be tested by comparing the gap closing with 
spin-spin correlations. Those can be numerically extracted 
by applying a strong magnetic field on both boundaries of 
the cylinder and measuring the AF order parameter (s^) 
as function of distance from the boundaries. This is shown 
in Fig. 4a for various values of Jh- For small Jh there is 
exponential decay as expected for the AKLT state [6] while 
for large Jh there is no decay as expected for the long 
range ordered Neel state. The extracted inverse AF spin- 
correlation length is plotted in Fig. 3b; its vanishing 
coincides with the vanishing of the gap, ruling out the 
possibility of an intermediate phase. The behaviour of the 
entanglement entropy Se when partitioning the system 
into two halves provides another signature of a single phase 


transition: a singularity develops at Jh = J, as seen in 
Fig. 4b. 

Boundary excitations- One of defining properties of 
the AKLT state are the emergent degrees of freedom on 
system boundaries. Consider a lattice vacancy obtained 
by removing a single site from the bulk of the honeycomb 
lattice; see Fig. 1(e). In the AKLT state, this breaks three 
valence bonds hence creating three unpaired spins-1/2 
surrounding the vacancy. Those are generically coupled, 
however, as long as the system respects time reversal sym¬ 
metry a Kramers-degenerate spin-1/2 degree of freedom 
must remain. This effect is reproduced in our numerical 
simulations of the model Eq. (3) with a vacancy: we 
remove a site and compute the ground state energy in 
sectors whose total differs by one. The light blue and 
purple curves in Fig. 3a represent the energy difference 
AA of those two states: they are degenerate in the AKLT 
phase. In contrast, the two states split in the magnetically 
ordered Neel phase, where AA becomes energy to create a 
magnon excitation. 

The free spin-1/2 degrees of freedom emerging at lattice 
vacancies in the AKLT phase can be detected experimen¬ 
tally due to their contribution to the free spin susceptibility, 
which will be proportional to the impurity density. 

Two site toy model - The mechanism by which inter¬ 
actions of the form of Eq. (1) arise from the directional 
nature of the orbitals may be intuitively understood by 
considering a two-site toy version of the model Eq. (3), 
H = -JH{J2-ySi,-f)^ - JH{Y.jS 2 ,-f)^ + Jsi^iS 2 ,i: see inset 
of Fig. 5. There is only one J coupling, which we take to 
be in the ^ = xy = \ orbital. 

The Hilbert space of three spins-1/2 at each site decom¬ 
poses into two doublets and a quartet. As we discussed 
before for large Jh/J the spin-3/2 quartets are energeti¬ 
cally favored by the Hund’s rule coupling, and they are 
coupled by an effective Heisenberg AF interaction Jeff = ^. 
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The Hilbert space of those two spins-3/2 decomposes in 
turn as follows: |(8)|=O©10203, and consequently 
the Heisenberg term may be written in terms of projectors 
to subspaces of total spin: 

= + (4) 

The energy splitting of those four multiplets is linear in 
J, as can be seen from dashed lines in bottom-left part of 
Fig. 5, where the spectrum of this 6 -spin model, computed 
via exact diagonalisation, is shown. As expected, the 
ground state is a singlet. Note, however, that this result is 
valid for small J/Jh only, when the doublets do not mix 
considerably with the spin-3/2 quartets. 

Upon increasing J/Jh the states belonging to doublets 
and quartets will mix; this can be seen in Fig. 5 in the form 
of level repulsion between the states of like total S'tot = 
0 , 1 , 2 (see color code) originating in | 0 |,| 0 | and | 0 | 
subspaces. Furthermore, the lowest S'tot = 0 , 1,2 levels 
are approximately degenerate - they form a groundstate 
manifold, while the unique state of Stot = 3 becomes 
an excited state. This mimics the spectral properties 
of the AKLT model Eq. (1) which on two sites reads 
F^aklt = in which all states but those forming a 

maximal total spin-3 between the two sites, are exactly 
degenerate. 



Figure 5. Eigenvalues of the two-site Hamiltonian as a function 
of J/Jh- Inset above: two-site toy model. 

The bunching of the lowest Stot = 0,1, 2 levels, a total of 
9 states, can be intuitively understood in the limit of large 
J/Jh, when a singlet forms between the spins coupled by 
J. The remaining two pairs of spins-1/2 form two spins-1 
due to the presence of Jh, whose Hilbert space is: 1 0 1 = 
0 0 10 2. The approximate degeneracy between these 
Stot = 0 , 1 , 2 states is broken once we include an isotropic, 
all orbital-to-all orbital coupling JJ2-y 7 ' ®i 7 ’ ^ 27 '- 

Adiabatic connection of band insulator and the AKLT 
state - We refer to the gapped phase in the diagram as 
the “AKLT” phase since at C/ = Jh = 0 the model can be 
smoothly connected in Hamiltonian space to the AKLT 


Hamiltonian, without closing the gap, as we now demon¬ 
strate. No other phase transitions intervening, the gapped 
phase of our diagram is therefore in the AKLT universality 
class. The result is similar in spirit to Refs. [12, 13] con¬ 
necting S = \ Heisenberg chain and the antiferromagnetic 
S' = 1/2 Heisenberg ladder, though in our intermediate 
transformations we break translation invariance. 

At the noninteracitng point the ground state of our 
model Eq. (2) is a product of decoupled singlets. We 
show that the AKLT state can be connected to such a 
product of singlets. The procedure works for AKLT spin 
models in arbitrary dimension. We work in the enlarged 
Hilbert space, where each AKLT spin-3/2 is composed of 
spin-1/2 degrees of freedom. The decoupling into singlets 
is performed repeating the following steps: (i) Choose a 
site i, and let {ij) for j = 1,2,3 denote the three links 
incident on i. (ii) Adiabatically weaken the coupling 
constants Jij , until site i is almost decoupled. The effective 
Hamiltonian for site % is given by spin-3/2 coupled to three 
effective spins-1/2 emerging at the boundary of the AKLT 
state around site i. (iii) This Hamiltonian is adiabatically 
connected to three decoupled singlets, each made up of 
one of the effective boundary spins- 1/2 and one of the 
three spins-1/2 which build the spin-3/2 on site i (see 
supplemental material for a generic example of a path in 
Hamiltonian space), (iv) Turn back the couplings. 
The three spins-1/2 at site i are now completely decoupled 
from each other (and only coupled to the spins-3/2 at 
neighbouring sites j). (v) Move to the next site. 

Upon visiting all sites, the Hamiltonian will have become 
a sum of decoupled singlets on individual links, without 
ever closing the gap. We have thus managed to find 
a path in hamiltonian space between Haklt and the 
noninteracting limit of our model Eq. (2). 

Conclusions and outlook - We considered states which 
can be realized in electronic systems with orbital degrees of 
freedom, in which electron hopping, due to the directional 
character of the orbitals, is highly anisotropic. For a half- 
filled 2D honeycomb lattice of t 2 g electrons we have shown 
that a 2D AKLT state may be realized as the ground 
state of the system. This state is known to be a universal 
resource for quantum computation [ 8 ]. 

The physics described in this paper, in particular the 
competition between the gapped singlet state of a valence 
bond solid (AKLT) type and that with a conventional long- 
range magnetic order, may be relevant for real correlated 
electron honeycomb systems. For example Na 2 lr 03 [14] 
and Na 2 Ru 03 [15] exhibit long-range magnetic order, but 
some samples of Li 2 Ru 03 have singlet valence bonds [16, 
17]. Though this state is at a different filling to the one 
we consider, it raises the hope to realize the VBS state in 
similar compounds. We note that we neglected the effective 
dd hopping via ligand (e.g. oxygen p-orbitals), which may 
be relevant in some materials; also the electron-lattice 
coupling may change the situation, leading to structural 
transitions with a loss of hexagonal symmetry, which seems 
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to be the case in M 0 CI 3 [18]. 

The general picture described above (possibility of for¬ 
mation of AKLT-like state in a multi-orbital system) can be 
realized also in other geometries, beyond the 2D honeycomb 
lattice: it may apply to recently synthesized hyperhoney¬ 
comb lattices or to a new 3D structure of 7 -Li 2 lr 03 [19], 
and to novel hyperoctagon lattice structures proposed by 
M. Hermanns and S. Trebst [20]. In all these cases, for 
ions in configuration, we anticipate the possibility of 
existence of two competing states, the gapless state with a 
long-range magnetic order, and the gapful state of a VBS 
type. 
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This work was supported by ISF and Marie Curie CIG 
grants (ES). 


[1] L. Balents, Nature 464, 199 (2010). 

[2] X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 
11413 (1989). 

[3] F. Haldane, Physics Letters A 93, 464 (1983). 

[4] F. D. M. Haldane, Phys. Rev. Lett. 50, 1153 (1983). 

[5] M. Hagiwara, K. Katsumata, 1. Affleck, B. 1. Halperin, 
and J. P. Renard, Phys. Rev. Lett. 65, 3181 (1990). 

[6] 1. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 
Rev. Lett. 59, 799 (1987). 


[7] 1. Affleck, T. Kennedy, E. Lieb, and H. Tasaki, Commu¬ 
nications in Mathematical Physics 115 , 477 (1988). 

[8] T.-C. Wei, 1. Affleck, and R. Raussendorf, Phys. Rev. Lett. 
106 , 070501 (2011). 

[9] J. c. V. Chaloupka, G. Jackeli, and G. Khaliullin, Phys. 
Rev. Lett. 105, 027204 (2010). 

[10] S. Lee, S. Choi, J. Kim, H. Sim, C. Won, S. Lee, S. A. 
Kim, N. Hur and J.-G. Park, J. Phys.: Condens. Matter 
24 , 456004 (2012). 

[11] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

[12] S. R. White, Phys. Rev. B 53, 52 (1996). 

[13] F. Anfuso and A. Rosch, Phys. Rev. B 75, 144420 (2007). 

[ 14 ] S. K. Choi, R. Coidea, A. N. Kolmogorov, T. Lancaster, 

I. 1. Mazin, S. J. Blundell, P. G. Radaelli, Y. Singh, 
P. Gegenwart, K. R. Choi, S.-W. Cheong, P. J. Baker, 
C. Stock, and J. Taylor, Phys. Rev. Lett. 108, 127204 
( 2012 ). 

[15] J. C. Wang, J. Terzic, T. F. Qi, F. Ye, S. J. Yuan, 
S. Aswartham, S. V. Streltsov, D. 1. Khomskii, R. K. 
Kaul, and G. Cao, Phys. Rev. B 90, 161110 (2014). 

[16] Y. Miura, Y. Yasui, M. Sato, N. Igawa, and K. Kakurai, 
Journal of the Physical Society of Japan 76 (2007). 

[17] G. Jackeli and D. Khomskii, Physical review letters 100 , 
147203 (2008). 

[18] H. Schafer, H.-G. V. Schnering, J. Tillack, F. Kuhnen, 
H. Wohrle, and H. Baumann, Z.Anorg.Allg.Chem. 353, 
281 (1967). 

[19] A. Biffin, R.D. Johnson, 1. Kimchi, R. Morris, A. Bombardi, 

J. G. Analytis, A. Vishwanath, R. Coidea, arXiv: 1407.3954 
(unpublished). 

[20] M. Hermanns, S. Trebst, Phys. Rev. B 89, 235102 (2014). 



arXiv: 1501.03512vl [cond-mat.str-el] 14 Jan 2015 


Supplemental material: Two-dimensional Valence Bond Solid (AKLT) states from t2g 

electrons 



Jh ( A = i ) 


fact that the lowest excitation carries spin-1 for both the 
AKLT and the Neel phase. 

Adiabatic connection While it is not essential, for read¬ 
ers’ interest, we show an explicit example of a path in 
Hamiltonian space mentioned in the main text, by means 
of exact diagonalization. The path itself is generic. 

We write the central spin-3/2 a.s S = 
the three effective boundary spins we denote by f^. The 
projector into the total spin of 3/2-1- 1/2 between S and 
one of the r spins is: 


Figure 1. Eigenvalues of H[Jh,A\ in Eq. (2) as function of Jh 
at A = 1, showing a singlet ground state for any Jh- 

DMRG simulations - We include here some tech¬ 
nical details about the DMRG simulations, performed 
using the developer branch of the ITensor package 
(http://itensor.org) on the Weizmann Institute cluster. We 
have extended the functionality of the package by writing 
a new matrix-product Hamiltonian file for a decorated 
honeycomb lattice with and without a vacancy (see Figs. 
ld,e in the main text). The computations were performed 
for cylindrical geometry, keeping the aspect ratio of cylin¬ 
der length to circumference (in units of hexagons) equal 
to 2. The system sizes simulated were 60, 126 and 216 
spin-1/2 degrees of freedom, which corresponds to 4x2, 
6x3 and 8x4 cylinders. The initial state for the DMRG 
algorithm was chosen to be a classical Neel state with 
^tot = 0- The free spin-1/2 degrees of freedom on the 
edges of the cylinder were frozen using a strong magnetic 
field applied there. These dangling free spins on the two 
edges of the cylinder live on two distinct sublattices of 
the Neel AF state, hence opposite magnetic fields were 
applied on different edges. We have kept up to m = 1600 
states in the 8x4 computation and up to m = 800 in the 
smaller ones. The gap was obtained by comparing the 
ground state energies in sectors with = 0,1 using the 


'P. + ^^{S,r,) = \{S + r,r-l. ( 1 ) 

Defining the interpolating Hamiltonian as: 

H[Jh, A] = -Jh ((E ^ 7 )^ - ^(^ + 1)) 

3 3 

+(1 ~ A) E '^§+5^'^’ ”^ 7 ) + E^^7 ■ W + (2) 

7=1 7=1 

the initial Hamiltonian is iLmit = R^[oo 7 0], and the final 
Hamiltonian is ffflnai = Ff[0,1]. A path in the (J//, A) 
parameter space should now connect these two points 
without closing the gap. 

A convenient path is via {Jh, A) = ( 00 ,0) —>■ { 00 , 1) —)■ 
(0,1). It is obvious that along the first transformation 
{Jh, A) = ( 00 , 0) —)■ ( 00 ,1) the ground state remains un¬ 
changed and the gap remains finite, since S is locked into a 
spin-3/2 state by Jh and the Hamiltonian only undergoes 
a change by a constant and by a positive proportional¬ 
ity coefficient. The results of the second transformation 
( 00 ,1) —>■ (0,1) is shown in Fig. 1, where the spectrum 
of the exactly diagonalized Hamiltonian is plotted as a 
function of Jh- the gap remains finite. We stress that this 
particular path was chosen for illustrative purposes only. 









